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INTRODUCTION 


Properly  designed  antenna  arrays  can  be  used  to  control  or  shape  radia¬ 
tion  patterns  for  covert  applications,  to  ensure  or  exclude  coverage  areas, 
and  to  reject  jamming  signals  and  reduce  other  electrical  interference.  The 
problem  of  array  design,  however,  is  particularly  difficult  when  the  antenna 
elements  are  closely  coupled  and  are  of  different  lengths.  Antenna  arrays  are 
usually  designed  by  trial  and  error  because  the  requirement  to  include  the 
effects  of  mutual  coupling  between  the  array  elements  complicates  design 
synthesis.  This  report  describes  an  approach  toward  design  synthesis  that 
combines  a  procedure  for  optimizing  feed  excitation  voltages  with  the  method 
of  moments  and  which  therefore  properly  includes  mutual  coupling. 

The  optimization  of  excitation  voltages  to  obtain  a  desired  radiation 
pattern  is  critical  to  the  design  of  many  antenna  systems,  especially  broad¬ 
cast  antenna  systems.  Increasingly  intense  competition  for  the  limited 
frequency  allocation  in  the  AM  broadcast  band  (550-1600  kHz)  has  led  to  the 
use  of  directional  array  antennas  to  reduce  interference.  In  this  band,  all 
available  channels  have  been  assigned.  For  a  new  station  to  receive  a 
frequency  assignment,  a  guarantee  must  be  made  that  there  will  be  no  inter¬ 
ference  with  existing  stations.  Thus,  new  stations  must  be  located  a 
significant  distance  away  from  any  existing  station  on  that  frequency.  Often 
additional  protection  must  be  provided  by  using  directional  array  antennas. 
In  the  United  States  more  than  25%  of  the  AM  broadcast  antennas  are  direc¬ 
tional  arrays.  This  allows  considerably  more  use  of  each  channel  and  is 
perhaps  the  first  attempt  at  practical  "frequency  reuse"  [1]. 

Typically,  the  directional  antenna  design  is  based  on  a  coverage  and 
protection  specification.  A  typical  coverage  specification  would  comprise  the 
following  conditions: 

(a)  Maximize  radiation  over  the  sector  covering  the  city  of  license. 

(b)  Maximize  radiation  at  headings  where  distant  locations  are  to  be 
protected. 

The  level  of  radiation  protection  (null  depth)  varies  with  frequency,  dis¬ 
tance,  and  the  power  level  of  the  local  transmitter  and  the  transmitter  to  be 
protected.  The  required  width  of  the  nulls  depends  upon  propagation  condi¬ 
tions  and  the  geometry  of  the  region  to  be  protected.  For  daytime  coverage, 
only  ground  wave  propagation  must  be  considered,  that  is,  only  nulls  at  an 
elevation  of  0  degrees.  For  the  night  time,  because  of  the  existence  of  sky 
wave  propagation,  elevation-pattern  nulls  must  be  considered.  The  directional 
array  antenna  design  must  necessarily  combine  both  quarter  wavelength  and 
taller  antenna  towers  to  achieve  elevation-pattern  nulls. 

Standard  practice  in  the  design  of  AM  broadcast  directional  array  an¬ 
tennas  has  been  to  calculate  the  array  excitation  voltages,  mutual  impedances, 
and  patterns  on  the  basis  of  an  assumed  sinusoidal  current  distribution  on 
each  antenna.  The  pattern  of  the  array  is  then  the  superposition  of  the 
fields  from  each  antenna  individually  (2J.  The  problem  with  this  approach  is 
that  the  currents  on  the  individual  antennas  are  not  necessarily  sinusoidal. 
Consequently,  the  actual  mutual  impedances  are  not  as  predicted.  The  actual 
field  for  each  antenna  (with  the  other  antennas  short-circui ted )  must  account 


for  scattering  from  the  other  antennas.  This  scattering  exists,  because  a 
current  distribution  can  exist  on  an  antenna  tower  even  though  the  input  is 
short-circuited.  The  sinusoidal  current  distribution  does  not  consider  such 
interactions.  In  general,  the  taller  the  AM  broadcast  antenna  towers,  the 

more  pronounced  will  be  the  difference  between  the  actual  array  patterns  and 

the  array  patterns  calculated  on  the  basis  of  sinusoidal  currents.  In  prac¬ 
tice,  the  sinusoidal  current  assumption  has  in  many  cases  led  to  incorrect 
designs  [3]. 

The  Mini -Numerical  Electromagnetics  Code  (MININEC )  1 4  J  is  a  computer 

program  which  calculates  accurate  antenna  characteristics,  including  actual 
current  distributions,  self-  and  mutual  impedances,  and  radiation  patterns. 
The  procedure  described  in  the  following  pages  for  the  optimization  of  the 
feed  excitation  voltages  of  an  antenna  system  to  yield  a  desired  radiation 

pattern  makes  use  of  MININEC  and  is  implemented  on  a  microcomputer. 

MININEC 

MININEC  is  a  frequency-domain  method-of -moments  computer  code  for  the 

analysis  of  thin-wire  antennas  [4].  MININEC  is  written  in  the  BASIC  language 
for  use  on  microcomputers  with  as  little  as  64K  of  memory.  The  MININEC 

program  is  based  on  the  numerical  solution  of  an  integral  representation  of 
the  electric  fields.  A  modified  Galerkin  procedure  is  used  to  solve  the 

integral  equation  [5],  This  formulation  results  in  a  uniquely  compact  code 
(that  is,  one  requiring  little  core  computer  storage)  suitable  for  implemen¬ 
tation  on  a  microcomputer.  The  MININEC  code  solves  for  the  impedance  and 

currents  on  arbitrarily  oriented  wires,  including  configurations  with  multiple 
wire  junctions.  Options  include  lumped  parameter  impedance  loading  of  the 
wires  and  the  calculation  of  radiation  patterns  in  free  space  and  over  flat, 
perfectly  conducting  ground  planes. 

Several  options  have  recently  been  added  to  MININEC  to  enhance  its 

utility  for  the  design  and  analysis  of  wire  antennas.  An  approximation  based 
on  the  Fresnel  plane-wave  reflection  coefficients  was  included  to  evaluate  the 
effect  of  imperfect  grounds  (impedance  not  equal  to  zero)  on  the  radiation 

pattern.  In  MININEC  the  current  is  computed  over  a  perfect  ground,  and  the 
radiated  field  can  then  be  computed  over  imperfect  grounds.  Multiple  media 
(that  is,  changes  in  surface  impedance)  can  be  specified.  The  choice  of 
medium  in  the  calculation  is  dependent  on  the  reflection  point.  A  radial 
ground  screen  approximation  (modified  reflection  coefficient)  can  also  be 

used.  Diffraction  from  the  edge  of  the  screen  is  not  included. 

The  analysis  of  an  antenna  problem  generated  by  a  thin-wire  method-of- 
moments  code  is  at  best  an  approximation.  Nonetheless,  once  the  inherent 

limitations  of  this  approach  are  taken  into  account,  highly  accurate  answers 
can  be  obtained  by  carefully  modeling  the  antenna  configuration. 

POWER  GAIN 

The  radiation  pattern  specifications  of  an  antenna  system  are  given  in 
terms  of  power  gain  or  directive  gain.  Both  specifications  are  given  in  terms 
of  decibels  (dB).  The  gain  for  the  direction  0  and  t>  is 


g  =  4if  P(e,<t>)/Pin  (i ) 

where  P(6,$)  is  the  power  radiated  per  unit  solid  angle  in  the  specified  di¬ 
rection.  P.  is  the  total  power  to  the  antenna  from  the  transmitter, 
in 

P.  =  4  R  (V»I*)  (2) 

in  2  e 

where  V  is  the  applied  source  voltage,  and  I*  is  the  conjugate  of  the  result¬ 
ing  feed  excitation  current.  In  addition 

2 

p(e,<i )  =  E*E*  (3) 

where  R  is  the  radius  of  the  observation  sphere;  E  is  the  far  electric  field; 
and  n  is  the  intrinsic  impedance  of  free  space. 

The  only  difference  between  directive  gain  and  power  gain  is  that  in  the 

latter  P  is  replaced  by  P  ,  where 
in  rad 

P  J  =  P.  -  P,  (4) 

rad  in  loss 

and  P,  is  power  loss  in  the  antenna  system.  The  gain  in  dB  is  expressed  as 
loss 

GdB  =  10  log  (G)  (5) 

For  the  purpose  of  feed  excitation  optimization,  the  multiple  feeds  of 
the  antenna  system  are  considered  to  be  an  N-port  network,  where  N  is  the 
number  of  feeds.  The  relationship  between  the  port  voltages  and  currents  can 
be  expressed  as 

N 

I .  =  y*  Y. .  V .  (6) 

i  n  i 

j=1 

where  Y,.  are  the  traditional  short-circuit  admittance  parameters.  The 
Y-parame\irs  are  defined  by  [6] 


Y. 


All  V  =  0 
k 

Except  k  =  j 


The  total  field  pattern  of  the  antenna  system  is  equal  to  the  sum  of  the 
fields  due  to  the  excitation  of  the  individual  feeds.  Changing  the  magnitude 
and  phase  of  the  individual  excitation  voltages  is  the  multiplication  of  the 
fields  from  each  individual  feed  (with  the  other  feeds  short-circuited)  by  the 
appropriate  amplitude  and  phase  prior  to  the  summing  for  the  total  pattern. 
Thus,  if  e  (9,  <p  )  is  the  far  electric  field  due  to  the  implied  voltages  of  V  * 
1,  the  total  far  electric  field  is 


3 


(8) 


E  ( 6,  «(> )  =  E  V.  e  (6,*) 
i  =  1  1  1 

where  V.  would  be  the  actual  feed  excitation  voltages  impressed  on  the  antenna 
system.1  Substituting  equations  (2),  (3),  (7),  and  (8)  into  equation  (1) 


G  ( 0 ,  $ )  - 


4  R  E  V  e  <e,*> 

i=1  1  L 


E  V  <e,*) 

i =1  1 


N  I  N  * 

“  R  I  V  ”  Y  y 


n  E  *e  1  E  ij  j 

i=1  L  3=1 


Thus  the  power  gain  can  be  calculated  for  any  set  of  feed  excitation 
voltages  e^  (6,  $)  and  the  short-circuit  admittance  parameters,  Y_,  are 

known.  These  values  can  be  calculated  using  MININEC. 

OPTIMIZATION  ALGORITHM 

For  the  purpose  of  optimizing  the  excitation  voltages  of  an  array  antenna 
to  provide  for  a  desired  radiation  pattern,  the  following  error  function  is 
defined: 

NP 

ERROR  =  E  W  (e,*)(R  (e  ,$  )-g.  (e  ,*  )  3  do) 

3  3  3 

where  NP  is  the  number  of  pattern  points  used  to  define  the  required  radiation 
pattern.  R^  (6,  )  is  the  required  radiation  pattern  in  dB  at  the  point  j. 

G^  (e,  #)  is  the  achieved  antenna  pattern,  which  is  given  by  equation  (9)  in 

dB.  The  error  is  defined  over  the  entire  synthesis  range.  The  least  squares 
error  is  used  to  ensure  that,  when  the  error  is  a  minimum,  each  component  of 
the  sum  is  a  minimum.  Therefore,  the  best  approximation  to  R_.  (6,  <t>  )  by  G_.(e, 

$  )  is  achieved  in  the  mean,  and  W_.(e,  $)  is  a  weighting  function  that  allows 

the  synthesis  precision  to  be  changed  over  the  synthesis  range. 

The  objective  of  optimization  is  to  minimize  the  error  function  of  equa¬ 
tion  (10).  A  BASIC  program,  called  OPTIMAL,  was  written  for  the  IBM  PC  to 
perform  this  optimization.  The  Appendix  contains  the  key  variables  and  a 
program  listing  of  OPTIMAL. 

OPTIMAL  uses  Newton's  method  to  perform  the  optimization  [7]  [8], 

Newton's  method  is  the  most  widely  used  of  all  iteration  formulas.  The  method 
uses  extrapolation  based  on  a  line  that  is  tangent  to  the  curve.  For  the  case 
of  a  single  variable,  the  method  can  be  developed  from  a  Taylor's  expansion  of 
the  form 


xn..  =  X  -  kF (X  )/F ’ (X  ) 

n+i  n  n  n 


4 


.vhere  X 


:  is  equivalent  to  the  point  where  the  curve  tangent 


through nthe  x  axis.  The  process  is  repeated  using 


x  =  x  , 

n  n+1 


(12) 


as  a  new  base  point.  When  the  value  of  F(X^+1  )  is  sufficiently  Small  the 
process  is  terminated.  The  progress  made  with  each  iteration  is  governed  by 


k,  the  gradient  constant.  If  k  is  chosen  to  be  too  small,  the  method  will 
converge  very  slowly.  If  k  is  too  large,  the  process  may  diverge.  In  OPTIMAL 
a  best  k  is  cnosen  using  a  type  of  binary  search  method  [7].  First,  the 
function  is  computed  for  k  equal  to  a  minimum  value  and  k  equal  to  a  maximum 
value.  In  the  case  of  OPTIMAL  the  initial  maximum  value  of  k  was  set  to 
.031  25,  and  the  initial  minimum  value  of  k  was  set  to  0.  A  new  function  is 
calculated  using  a  value  of  k  which  is  halfway  between  the  two  k  values.  If 
the  new  function  is  greater  than  the  function  with  the  maximum  k  value  and  the 
function  with  the  maximum  k  is  less  than  the  function  with  the  minimum  k 
value,  the  minimum  k  is  changed  to  the  mean  or  halfway  k.  If  the  new  function 
is  less  than  the  function  with  the  maximum  k  value  and  the  maximum  k  function 
is  greater  than  the  function  with  the  minimum  k  value,  the  maximum  k  value  is 
changed  to  the  mean  or  halfway  k  value.  The  process  is  repeated  until  the 
interval  between  the  minimum  and  maximum  k  values  reaches  some  minimum  value. 
In  the  case  of  OPTIMAL  this  minimum  value  is  .0000001. 


The  error  function  of  equation  (10)  is  a  function  of  2N  parameters. 
There  is  both  a  real  and  an  imaginary  voltage  for  each  of  the  feed  points. 
Using  V  to  denote  each  of  these  2N  parameters,  each  successive  iteration  of 
all  2N  parameters  is  given  by 


Error  (V  ) 

V  =  V - — 

n+1  n  d [ Error ( V  )] 


(13) 


n 


dV 


Analytical  expressions  are  used  for  each  of  the  required  derivatives  in 
OPTIMAL. 


SAMPLE  PROBLEM 


MININEC  was  used  to  model  an  array  of  three  antenna  elements.  The  short- 
circuit  admittance  parameters  were  thus  calculated  for  feeds  on  eacn  of  tne 
three  antenna  elements.  For  each  calculation,  a  single  feed  was  impressed 
with  a  voltage  of  1+]j<)  with  the  other  two  feeds  short-circuited.  The  far 
electric  field  was  then  calculated  tor  each  of  the  three  feeds.  Four  pattern 
points  for  each  of  tne  antenna  feeds  were  computed  for  this  example.  Addi¬ 
tional  or  fewer  pattern  points  may  be  chosen  depending  on  the  application.  As 
described  above,  these  are  the  required  inputs  to  tne  program  OPTIMAL.  OPTIMAL 
will  read  these  data  from  an  input  file  and  write  the  results  to  an  output 
tile.  The  following  is  the  input  file  that  was  used  for  this  sample  problem. 
Comments  have  been  added  in  italics  to  identify  the  data. 


3 


■.Number  of  feed  points 
•.Short  circuit  parameters 


. 1 20627E-1 ,-.1771 22E-1 
.241 459E-2, - • 1 64955E-1 
•639456E-3 , - . 1 4609E-2 
.261 697E-2, - . 256524E-1 
.241 457E-2, . 1 64956E-1 
. 1 20627E-1 ,-.1771 22E-1 
4 
0 

2.158E-1 ,-.380 
-.893E-1 ,1 .239 
-.557, 1.428 
-.981 , 1 .446 
6.233E-1 , 3.055E-1 
-.246,4.334E-1 
1 . 491 E-1 , 3 . 752E-1 
6.233E-1 , 3.055E-1 
-.981 , 1 .446 
3.952E-1 ,4.996E-1 
3.758E-1 ,5.21 3E-3 
2.1 58E-1 ,-.380 
10,1 
-11,1 
-40,  1 
-35,1 
250 
4 
Y 
10 


Number  of  pattern  points 
Radius  of  observation  sphere 
Far  electric  fields 


■Required  gains  and  weights 


Maximum  number  of  iterations 
Minimum  least  squares  error 
Option  for  listing  errors 
Number  of  iterations  before  listing 


The  output  file  which  results  from  the  above  input  file  is  given  below. 
Editorial  remarks  that  have  been  added  are  in  italics. 


OPTIMAL 

IBM  PC  VERSION 
10-24-1985 


The  maximum  number  of  feed  points  is  10.  This  number  can  be  increased  bg 
changing  the  arrays  of  lines  120  -  130  of  OPTIMAL  to  the  appropriate  value. 
In  addition  the  flag  of  line  6b0  would  have  to  be  changed. 


NUMBER  OF  FEED  POINTS?  3 
ANTENNA  ADMITTANCE  VALUES 


REAL 

AND 

IMAG 

ADMITTANCE 

Y( 

1 

/ 

i 

)? 

.0120627  ,-. 

01 77122 

REAL 

AND 

IMAG 

ADMITTANCE 

Y  ( 

1 

t 

2 

)? 

2.41 459E-03 

,  .0164955 

REAL 

AND 

IMAG 

ADMITTANCE 

Y  ( 

1 

r 

3 

)? 

6. 39456E-04 

,-.001 4609 

REAL 

AND 

IMAG 

ADMITTANCE 

Y  ( 

2 

t 

2 

)? 

2.61 697E-03 

,  .0256524 

REAL 

AND 

IMAG 

ADMITTANCE 

Y  ( 

2 

f 

3 

)? 

2.41 457E-03 

,  .0164956 

REAL 

AND 

IMAG 

ADMITTANCE 

Y  ( 

3 

t 

3 

)? 

.0120627 

01 771 22 

The  maximum  number  of  pattern  points  is  20.  This  number  ean  be  increased 
by  changing  the  arrays  of  lines  120-130  of  OPTIMAL  to  the  appropriate  value. 
In  addition,  the  flag  of  line  920  would  have  to  be  changed. 


NUMBER  OF  PATTERN  POINTS?  4 

When  the  observation  sphere  radius  is  set  equal  to  zero ,  the  effect  of 
range  is  ignored.  This  is  similar  to  the  convention  used  in  MININEC  in  the 
calculation  of  pattern  gains. 

OBSERVATION  SPHERE  RADIUS  (M)?  0 

PATTERNS  ( REAL, IMAG)  OF  INDIVIDUAL  FEEDS  (V/M) 

IT  IS  ASSUMED  THAT  FEED  VOLTAGES  WERE  (1,0). 


PATTERN 

POINT 

1 

OF 

FEED 

POINT 

1 

7 

.2158  ,-.38 

PATTERN 

POINT 

2 

OF 

FEED 

POINT 

1 

7 

-.0893  ,  1.239 

PATTERN 

POINT 

3 

OF 

FEED 

POINT 

1 

7 

-.557  ,  1.428 

PATTERN 

POINT 

4 

OF 

FEED 

POINT 

1 

7 

-.981  ,  1.44b 

PATTERN 

POINT 

1 

OF 

FEED 

POINT 

2 

7 

.6233  ,  .3055 

PATTERN 

POINT 

2 

OF 

FEED 

POINT 

2 

7 

-.246  ,  .4334 

PATTERN 

POINT 

3 

OF 

FEED 

POINT 

2 

7 

.1491  ,  .3752 

PATTERN 

POINT 

4 

OF 

FEED 

POINT 

2 

7 

.6233  ,  .3055 

PATTERN 

POINT 

1 

OF 

FEED 

POINT 

3 

7 

-.981  ,  1.446 

PATTERN 

POINT 

2 

OF 

FEED 

POINT 

3 

7 

.3952  ,  .4996 

PATTERN 

POINT 

3 

OP 

FEED 

POINT 

3 

7 

.3758  ,  .005213 

PATTERN 

POINT 

4 

OF 

FEED 

POINT 

3 

7 

.2158  ,-.38 

DESIRED  RADIATION  PATTERN  (DB) 

PATTERN  AND  WEIGHT  OF  POINT  1  ?  10  ,  1 

PATTERN  AND  WEIGHT  OF  POINT  2  ?  -11  ,  1 

PATTERN  AND  WEIGHT  OF  POINT  3  ?  -40  ,  1 

PATTERN  AND  WEIGHT  OF  POINT  4  ?  -35  ,  1 

NUMBER  OF  ITERATIONS?  25( 

MINIMUM  LEAST  SQUARES  ERROR?  4 

LISTING  OF  ITERATIONS  AND  ERROR  (Y/N)?  Y 
NUMBER  OF  ITERATIONS  BEFORE  LISTING?  10 


Initially,  jPTI'-'AL  is  slow  in  converging.  However,  once  it  starts  to 
ion  verge ,  it  converges  prickly.  On  the  IiiM  PC  the  time  of  the  following  cal¬ 
culation  is  greater  than  lb  minutes.  The  time  of  any  calculation  is  dependent 
on  the  number  >f  fee.:  points  and  pattern  points.  The  number  of  iterations 
require a  for  convergence  is  legenaent  on  the  word  length  of  the  calculations. 
For  shorter  won:  lengths,  it  wi . i  take  more  iterations  for  convergence. 


m 


1 


ITERATION 

ERROR 

1  0 

1 233.41 1 

2  1 

1 222.322 

is 

1219.684 

•*  \ "  • 

40 

1197.274 

■>  i 

1  193.403 

'■*  *.'•»  *  * 

hO 

1 1 88.957 

■>v-; 

/u 

1 184.84 

m 

BO 

1  182.973 

JO 

1181.16 

7 


100 

1178.017 

110 

1172.157 

120 

1165.476 

1  30 

1 1 59.169 

140 

1157.159 

150 

1154.309 

160 

1 1 51 .327 

170 

1149.454 

180 

1142.514 

190 

1 1 34.836 

200 

887.8659 

210 

382.5797 

220 

3.573898 

**********************  OUTPUT  ************************ 

The  minimum  least  squares  error  is  achieved  in  220  trials. 

NUMBER  OF  ITERATIONS  =  220 


LEAST  SQUARES  ERROR  =  3.573898 


*****  SYNTHESIZED  FEED  POINT  VOLTAGES  ***** 

FEED  NO.  REAL  (V)  IMA G  (V)  MAGNITUDE  (V)  PHASE  (DEG) 

1  -7.468426E-02  .2037476  .2170042  110.1306 

2  1.093298  1.2362  1.6503  48.51037 

3  1.480106  -1.013934  1.794095  -34.41275 

The  above-synthesized  feed  point  voltages  were  used  in  MININEC  and 
confirmed  the  radiation  pattern  predicted  below. 

*****  RADIATION  PATTERN  (DB)  ***** 

PATTERN  NO.  REQUIRED  (DB)  SYNTHESIZED  (DB) 

1  10  11.08301 

2  -11  -9.60534 

3  -40  -40.50869 

4  -35  -34.55599 

INPUT  POWER  (WATTS)  =  2.435024E-02 

It  is  possible  with  OPTIMAL  to  continue  the  synthesis  process. 

SUMMARY  AND  RECOMMENDATIONS 


A  microcomputer  can  use  the  program  MININEC  and  OPTIMAL  to  optimize  the 
antenna  system  feeds  for  a  desired  radiation  pattern.  This  capability 
demonstrates  the  growing  importance  of  the  microcomputer  in  antenna  modeling. 

Several  areas  for  further  study  can  be  suggested.  The  Newton  method  is 
slow  and  does  not  guarantee  that  local  minimums  are  being  avoided.  More 
modern  methods  of  optimization  should  be  considered  for  this  problem.  A 
conjugate  gradient  method  such  as  described  by  Fletcher  and  Powell  [9]  or  the 
new  least  squares  algorithm  of  Dennis  et  al.  [10]  are  techniques  that  should 
be  considered.  In  this  study  only  optimization  of  feed  excitation  voltages 
was  considered.  It  is  also  possible  to  optimize  antenna  locations 
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and  heights.  Such  an  optimization  capability  would  approach  a  total  design 
capability  for  the  broadcast  directional  array  antenna  system  and  may  have 
useful  applications  to  other  directional  antenna  array  systems. 
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APPENDIX 


I 

!  KEY  VARIABLES  AND  PROGRAM  LISTING  OF  OPTIMAL 


Included  in  this  Appendix  are  the  definitions  of  the  key  variables  and 
the  program  listing  of  the  BASIC  program  OPTIMAL. 

LISTING  OF  KEY  VARIABLES 


A  (I ) 

AI ( I , J  ) 

AR  (I,J) 

CG 

C$ 

SKI,  J) 

ER(I, J ) 

FI  (I  ) 
FR(I) 

F$ 

GC 

G  ( I ,  1  ) 
G(I, 2 ) 
IT 


Achieved  gain  (dB)  at  pattern  point  I 

Imaginary  value  of  far  electric  fields  (V/m )  from  feed  point  I 
with  value  of  (1,0)  at  pattern  point  J 

Real  value  of  far  electric  field  (V/m)  from  feed  point  I  with 
value  of  (1,0)  at  pattern  point  J 

Maximum  constant  for  gradient 

Flag  for  continuation  of  iterations 

Imaginary  value  of  far  electric  field  (V/m)  from  feed  point  I  for 
actual  voltages  at  pattern  point  J 

Real  value  of  far  electric  field  (V/m)  from  feed  point  I  for 
actual  voltages  at  pattern  point  J 

Imaginary  value  of  voltage  at  feed  point  I  before  iterations 
Real  value  of  voltage  at  feed  point  I  before  iteration 
Names  of  files 

(2*  pi  )/ (intrinsic  free  space  impedance) 

Least  squares  derivative  for  real  value  of  feed  point  I 
Least  squares  derivative  for  imaginary  value  of  feed  point  I 
Number  of  iterations 


KO 

1$ 

L 

LR 

MG 

NF 


Minimum  interval  for  binary  search 

Flag  for  using  input  file 

Changes  natural  log  to  log  base  10 

Log  to  base  10  of  sphere  radius  squared  (dB) 

Magnitude  of  feed  point  voltage  (V) 

Number  of  feed  points 


NI 

NL 

NP 

0$ 

P 

POWER 
PO 
RAD 
R  (I  ) 

SO 

T(I) 

TT  ( I  ) 

VI  (I) 
VR(I) 
W(I  ) 

*KI, I) 

YR  ( I , I  ) 


Number  of  iterations 

Number  of  iterations  between  listing  errors 

Number  of  pattern  points 

Flag  for  using  output  file 

Phase  of  feed  excitation  voltage  (deg) 

Input  power  (W ) 

Changes  degrees  to  radians 
Observation  sphere  radius  (m) 

Required  gain  (dB)  at  pattern  point  I 
Least  squares  error 

Difference  in  required  and  achieved  gain  at  pattern  point  I  (dB ) 

Magnitude  of  far  electric  field  (V/m)  due  to  actual  feed 
excitation  voltages  at  pattern  point  I 

Imaginary  value  of  synthesized  voltage  at  feed  point  I 

Real  value  of  synthesized  voltage  at  feed  point  I 

Weighted  value  of  least  squares  error  at  pattern  point  I 

Imaginary  part  of  short-circuit  admittance  (mhos) 

Real  part  of  short-circuit  admittance  (mhos) 


PROGRAM  LISTING 


The  following  BASIC  program,  OPTIMAL,  was  written  to  run  on  an  IBM-PC. 
However,  OPTIMAL  has  been  run  on  several  IBM-compatible  computers.  With  only 
minor  changes  it  should  be  possible  to  run  OPTIMAL  on  most  computers. 


100  REM  -  NUMBER  OF  FEED  POINTS  =  10 

110  REM  -  NUMBER  OF  PATTERN  POINTS  =  20 

120  DIM  A( 20) , AR( 10, 20) , AI ( 10, 20) , ER( 20) , El ( 20) ,FR( 10) ,FI ( 10) ,G(10, 2) 
130  DIM  R( 20 ) , W( 20 ) ,T(20) , TT (20) ,VR(10) ,VI(10) ,YR(10, 10) ,YI(10, 10) 

140  REM  **********  DATA  INITIALIZATION  ********** 

150  REM  -  NO  CONTINUATION  OF  ITERATIONS 

160  C$="N" 

170  REM  -  CHANGES  RADIANS  TO  DEGREES 

180  P0=ATN(l)/45 

190  REM  -  CHANGES  NATURAL  LOG  TO  LOG  BASE  10 

200  L= 10/ LOG ( 10 ) 

210  REM  -  MAXIMUM  CONSTANT  FOR  GRADIENT 

220  CG=. 03125 

230  REM  -  MINIMUM  INTERVAL  FOR  BINARY  SEARCH 

240  K0=. 0000001 

250  REM  -  (2  *  PI )/( INTRINSIC  FREE  SPACE  IMPEDANCE) 

260  GC=. 016678 

270  REM  **********  INPUT/OUTPUT  FILES  ********** 

280  PRINT  "INPUT  FILE  (Y/N) 

290  INPUT  1$ 

300  IF  I $="N"  THEN  350 
310  IF  I  $  <  > " Y"  THEN  280 
320  PRINT  "NAME  OF  FILE 
330  INPUT  F$ 

340  OPEN  F$  FOR  INPUT  AS  *  1 
350  PRINT  "OUTPUT  FILE  (Y/N)"; 

360  INPUT  OS 

370  IF  0$= "N"  THEN  420 

380  IF  0$  <  > "Y"  THEN  350 

390  PRINT  "NAME  OF  FILE  "; 

400  INPUT  F$ 

410  OPEN  F$  FOR  OUTPUT  AS  #2 

420  REM  **********  TITLE  ********** 

430  W$  =  "  *****************- 

440  X$="  OPTIMAL" 

450  Y$  = "  IBM  PC  VERSION" 

460  Z$="  " +DATES 


470 

PRINT 

ws 

400 

PRINT 

X$ 

490 

PRINT 

YS 

500 

PRINT 

Z$ 

5  10 

PRINT 

WS 

520 

PRINT 

5  30 

IF  OS 

=  "  Y  " 

THEN 

PRINT 

a  2 ,  w$ 

540 

IF  OS 

=  "  Y" 

THEN 

PRINT 

#2  ,  XS 

550 

IF  OS 

= »»  y  " 

THEN 

PRINT 

S2  ,  YS 

560 

IF  OS 

_  H  Y» 

THEN 

PRINT 

#2  ,  ZS 

1  2 


A  Al\A 


570  IF  0$  =  "Y"  THEN  PRINT  #2,W$ 

580  IF  0§=" Y"  THEN  PRINT  #2, 

590  REM  **********  DATA  INPUT  ********** 

600  Q$="?  " 

610  REM  -  ANTENNA  SHORT-CIRCUIT  ADMITTANCE  VALUES 

620  W$  =  "  NUMBER  OF  FEED  POINTS'* 

630  PRINT  W$ ; 

640  IF  I $= "N"  THEN  INPUT  NF  ELSE  INPUT  #1,NF 
650  IF  NF< 11  THEN  690 
660  PRINT 

670  PRINT  "NUMBER  OF  FEED  POINTS  EXCEEDS  DIMENSION..., 
680  IF  I $=" Y"  THEN  STOP  ELSE  620 
690  IF  0$="Y"  THEN  PRINT  #2,W$+Q$?NF 
700  W$="  ANTENNA  ADMITTANCE  VALUES" 

710  PRINT  W$ 

720  IF  0$= "Y"  THEN  PRINT  #2,W$ 

730  FOR  1=1  TO  NF 
740  FOR  J=I  TO  NF 

750  W$= "REAL  AND  IMAG  ADMITTANCE  Y( " 

760  PRINT  W$ ; I ; “ , " ; J ; " ) " ; 

770  IF  I$="N"  THEN  INPUT  YR { I , J ) , YI ( I , J ) 

780  IF  I $= " Y"  THEN  INPUT  # 1 , YR ( I , J ) , YI ( I , J ) 


790  IF  0$="Y"  THEN  PRINT  #2 , W$ ; I ; " , " ; J ; " ) " ?Q$ j YR( I , , 

800  REM  -  THE  ADMITTANCE  VALUES  ARE  SYMMETRICAL 

810  IF  I=J  THEN  840 
820  YR ( J , I ) =YR ( I , J ) 

830  YI ( J , I ) =YI ( I , J ) 

840  NEXT  J 
850  NEXT  I 
860  PRINT 

870  IF  0$="Y"  THEN  PRINT  #2, 

880  REM  -  FAR  ELECTRIC  FIELD  FOR  UNIT  VOLTAGES 


)  " ; Q$  * YR( I,J);"," ?YI(I,J) 


890  W$=" 


NUMBER  OF  PATTERN  POINTS" 


900  PRINT  W$; 

910  IF  I$="N"  THEN  INPUT  NP  ELSE  INPUT  #1,NP 
920  IF  NP<21  THEN  960 
930  PRINT 

940  PRINT  "NUMBER  OF  PATTERN  POINTS  EXCEEDS  DIMENSIONS..." 

950  IF  I $=" Y"  THEN  STOP  ELSE  910 
960  IF  0$="Y"  THEN  PRINT  #2,W$+Q$;NP 

970  REM  -  OBSERVATION  SPHERE  RADIUS 

980  W$  =  "OBSERVATION  SPHERE  RADIUS  (M)" 

990  PRINT  W$ ; 

1000  IF  I $  = " N "  THEN  INPUT  RAD  ELSE  INPUT  tl , RAD 
1010  IF  0$="Y"  THEN  PRINT  #2,W$+Q$;RAD 

1020  REM  -  LOG  VALUE  TO  BE  USED  IN  ACHIEVED  GAIN  CALCULATION 

1030  LR=0 

1040  IF  RAD> 0  THEN  LR=2*L*LOG( RAD) 

1050  W$  =  "  PATTERNS  (REAL,  IMAG)  OF  INDIVIDUAL  FEEDS  (V/M)" 

1060  PRINT  W$ 

1070  IF  0$  = " Y"  THEN  PRINT  #2,W$ 

1080  W$  =  "  IT  IS  ASSUMED  THAT  FEED  VOLTAGES  VJERE  (1,0)." 


1000  IF  1$: 
1010  IF  0$- 
1020  REM  -■ 
1030  LR=0 
1040  IF  RA1 
1050  W$  =  " 
1060  PRINT 


1070  IF  0$= " Y" 
1080  W$="  IT 
1090  PRINT  W$ 
1100  IF  0$  =  " Y " 


THEN  PRINT  #2,W$ 


»*  */  ^  ■  u.*.  •  , 


V 


mo 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 


FOR  1=1  TO  NF 
FOR  J=1  TO  NP 
W$=" PATTERN  POINT  " 

X$  =  "  OF  FEED  POINT  " 

PRINT  W$; J;X$; I; 

IF  I$="N"  THEN  INPUT  AR( I , J ) , AI ( I , J ) 

IF  I$="Y"  THEN  INPUT  # 1 , AR ( I , J ) , AI ( I , J ) 

IF  0$  =  "Y"  THEN  PRINT  #2 , W$ j J ; X$ ; I ; Q$ ; AR( I , J ) ; “ , " ? AI ( I , J ) 

NEXT  J 
NEXT  I 

REM  -  DESIRED  RADIATION  PATTERN  ( POWER  GAIN) 

PRINT 

IF  0$= "Y"  THEN  PRINT  #2, 

W$="  DESIRED  RADIATION  PATTERN  (DB)" 

PRINT  W$ 

IF  0$ =  "  Y"  THEN  PRINT  #2,W$ 

FOR  1=1  TO  HP 

W$  = " PATTERN  AND  WEIGHT  OF  POINT  " 

PRINT  W$ ; I ? 

IF  I $  = " N "  THEN  INPUT  R(I),W(I)  ELSE  INPUT  #1,R(I),W(I) 

IF  0$  =  "Y"  THEN  PRINT  # 2 , W$ ; I ; Q $ ; R ( I ) ? " , " ; W ( I ) 

NEXT  I 
PRINT 

IF  0$  = " Y"  THEN  PRINT  #2, 

REM  -  NUMBER  OF  ITERATIONS 

W$="  NUMBER  OF  ITERATIONS" 

PRINT  W$; 

IF  I$  =  "N"  THEN  INPUT  NI  ELSE  INPUT  #1,111 
IF  0$=" Y"  THEN  PRINT  #2,W$+Q$;NI 
IF  C$="Y"  THEN  NI=NI+IT 

REM  -  MINIMUM  LEAST  SQUARES  ERROR 

W$="  MINIMUM  LEAST  SQUARES  ERROR" 

PRINT  W$; 

IF  I  $  =  "N"  THEN  INPUT  ER  ELSE  INPUT  #1,ER 
IF  0$="Y"  THEN  PRINT  #2,W$+Q$;ER 

REM  -  OPTION  OF  LISTING  ITERATIONS  AND  ERROR 

W$= " LISTING  OF  ITERATIONS  AND  ERROR  (Y/N)" 

PRINT  W$ ; 

IF  I$="N"  THEN  INPUT  L$  ELSE  INPUT  #1,L$ 

IF  0$=" Y"  THEN  PRINT  #2,W$+Q$;L$ 

IF  L$="N"  THEN  1570 

W$="  NUMBER  OF  ITERATIONS  BEFORE  LISTING" 

PRINT  W$; 

IF  I$="N"  THEN  INPUT  NL  ELSE  INPUT  #1,NL 
IF  0$  = " Y"  THEN  PRINT  #2,W$+Q$;NL 
IF  NL=0  THEN  L$="N" 

PRINT 

IF  0$  =  "Y"  THEN  PRINT  #2, 

IF  L$="Y"  THEN  PRINT  "  "," ITERATION" , "ERROR" 

IF  L$  =  "Y"  AND  0$  =  "Y"  THEN  PRINT  #2,"  "," ITERATION" ERROR" 

IF  C$="Y"  THEN  1780 

REM  *********  BEGIN  ITERATIONS  ********** 

IT=0 

REM  -  IF  FEED  VOLTAGES  ARE  NOT  (1,0)  FOR  PATTERN,  CHANGE 


ITERATION" , "ERROR" 

PRINT  #2,"  ", "ITERATION" , "ERROR" 


********** 


*v.v.v.v. 


1650  FOR  1=1  TO  NF 
1660  VR( I)=l 
1670  VI(I)=0 
1680  NEXT  I 

1690  REM  -  COMPUTE  FIRST  ACHIEVED  GAIN  AND  ERROR 

1700  GOSUB  3060 

1710  REM - FIRST  ERROR 

1720  S0=S 

1730  REM  **********  OUTPUT  CHECK  ********** 

1740  REM  -  IF  ITERATIONS  EQUAL  MAXIMUM, 

1750  IF  IT=NI  THEN  2380 

1760  REM  -  IF  COMPUTED  ERROR  IS  LESS  THAN  MINIMUM  ERROR, 

1770  IF  SO <  =ER  THEN  2380 

1780  REM  **********  ERROR  DERIVATIVE  ********** 

1790  REM  -  CONSIDER  ALL  FEED  POINTS 

1800  FOR  K=1  TO  NF 

1810  REM  -  CONSIDER  BOTH  REAL  AND  IMAGINARY  VOLTAGES 

1820  FOR  11=1  TO  2 
1830  D1 =0 

1840  REM  -  INPUT  POWER  DERIVATIVE 

1850  FOR  J=1  TO  NF 

1860  REM  -  FUNCTION  OF  REAL  VOLTAGE 

1870  IF  11=1  THEN  D1=D1 +VR ( J ) * YR ( K , J ) 

1880  REM  -  FUNCTION  OF  IMAGINARY  VOLTAGE 

1890  IF  11=2  THEN  D1 =D1 +VI ( J ) * YR ( K , J ) 

1900  NEXT  J 
1910  D=0 

1920  FOR  1=1  TO  NP 

1930  REM  -  ELECTRIC  FIELD  DERIVATIVES 

1940  IF  11=1  THEN  D2=2* ( ER{ I ) *AR( K, I ) +EI ( I ) *AI ( K, I ) ) 

1950  IF  11=2  THEN  D2=2* (El ( I ) *AR(K, I )-ER( I ) *AI ( K, I ) ) 

1960  D=D+W( I ) *T( I )*(Dl/POWER-D2/TT( I ) ) 

1970  NEXT  I 

1980  REM  -  LEAST  SQUARES  DERIVATIVE 

1990  G ( K , 1 1 ) =2*L*D 
2000  NEXT  II 
2010  NEXT  K 

2020  REM  *********  NEWTON  METHOD  ********** 

2030  REM  -  SAVE  PRESENT  VOLTAGES 

2040  FOR  1=1  TO  NF 
2050  FR ( I ) =VR ( I ) 

2060  F I ( I ) =  V I (  I  ) 

2070  NEXT  I 


2080  REM  -  ITERATIVE  METHOD  OF  FINDING  GAIN  CONSTANTS 

2090  C0=0 
2100  C=CG 

2110  REM  -  ERROR  CALCULATION  FOR  FIRST  GRADIENT  CONSTANT 

2120  GOSUB  3010 
2130  C2=CG 
2140  S2=S 

2150  REM - R I  NARY  SEARCH  FOR  REST  GRADIENT  CONSTANT 

2160  C=(C2+C0)/2 

2170  REM  -  ERROR  CALCULATION  FOR  SUCCESSIVE  GRADIENT  CONSTANTS 

2180  GOSUB  3010 


2190  S1=S 
2200  C1=C 

2210  IF  S2 < SI  AND  S0>S2  THEN  2250 

2220  C2=C1 

2230  S2=S1 

2240  GOTO  2270 

2250  C0=C1 

2260  S0=S1 

2270  IF  C2-C0>K0  THEN  2160 


2280  REM  -  ERROR  CALCULATION  AFTER  BINARY  SEARCH  FOR  GRADIENT  CONSTANT 

2290  S0=S1 

2300  REM  *********  UPDATE  ITERATION  ********** 

2310  IT= IT+1 

2320  REM  -  ITERATION  AND  ERROR  LISTING  OPTION 


2330  IF  L$="N"  THEN  1750 
2340  TI=NL*INT( IT/NL) 

2350  IF  IT=TI  THEN  PRINT  "  " ,IT,S0 

2360  IF  IT=TI  AND  0$="Y"  THEN  PRINT  #2,"  ",IT,S0 

2370  GOTO  1750 

2380  REM  ********  OUTPUT  ********** 

2390  *******************  OUTPUT  ********************* 

2400  PRINT  W$ 

2410  PRINT 

2420  IF  OS  =  "Y"  THEN  PRINT  #2,VJ$ 

2430  IF  0$  = " Y"  THEN  PRINT  #2, 

2440  REM  -  NUMBER  OF  ITERATIONS 

2450  W$=" NUMBER  OF  ITERATIONS  =" 

2460  PRINT  W$ ; IT 

2470  IF  0$  = " Y"  THEN  PRINT  #2,W$;IT 

2480  REM  -  LEAST  SQUARES  ERROR 

2490  W$  =  " LEAST  SQUARES  ERROR  =" 

2500  PRINT  W$;S0 

2510  IF  0$="Y"  THEN  PRINT  #2,W$?S0 
2520  PRINT 

2530  IF  0$= " Y"  THEN  PRINT  #2, 

2540  REM  -  SYNTHESIZED  FEED  POINT  VOTAGES 

2550  W$  =  "***★*  SYNTHESISED  FEED  POINT  VOLTAGES  *****" 

2560  PRINT  "  " ,  W$ 

2570  IF  05  = " Y"  THEN  PRINT  #2,"  ",W$ 

2580  U$  =  " FEED  NO." 

2590  W$  =  " REAL  (V) " 

2600  X$=" IMAG  (V) " 

2610  Y$= "MAGNITUDE  (V)” 

2620  Z$= " PHASE  (DEG)" 

2630  PRINT  U$,W$,X$,Y$, Z$ 

2640  IF  0$=" Y"  THEN  PRINT  # 2 , U$ , W$ , X$ , Y$ , Z $ 

2650  FOR  1=1  TO  NF 

2660  REM  -  MAGNITUDE  AND  PHASE  CALCULATIONS 

2670  MG=SQR ( VR ( I ) *VR( I ) +VI ( I ) *VI ( I ) ) 

2680  IF  VR ( I ) <  >  0  THEN  2710 

2690  P=0 

2700  GOTO  2730 

2710  P=ATN (VI(I)/VR(I))/P0 


2720  IF  VR( I ) <0  THEN  P=P+SGN( VI ( I ) ) *180 
2730  PRINT  I, VR< I ) , VI ( I ) ,MG, P 

2740  IF  0$= "Y"  THEN  PRINT  #2 , I , VR( I ) , VI ( I ) , MG, P 
2750  NEXT  I 
2760  PRINT 

2770  IF  0$  = "Y"  THEN  PRINT  #2, 


2780  REM  -  RADIATION  PATTERNS 

2790  W$="  *****  RADIATION  PATTERN  (DB)  *****" 


2800  PRINT  W$ 

2810  IF  0$  =  "Y"  THEN  PRINT  #2,W$ 

2820  W$=" PATTERN  NO." 

2830  X$  =  " REQU I RED  (DB)" 

2840  Y$=" SYNTHES I ZED  (DB)" 

2850  PRINT  W$,X$,Y$ 

2860  IF  0$=" Y"  THEN  PRINT  #2,W$,X$,Y$ 

2870  FOR  1=1  TO  NP 
2880  PRINT  I , R( I ) , A( I ) 

2890  IF  0$= "Y"  THEN  PRINT  #2 , I , R( I ) , A{ I ) 

2900  NEXT  I 

2910  REM - INPUT  POWER 

2920  W$  =  "  INPUT  POWER  (WATTS) =" 

2930  PRINT  W$  •  POWER 

2940  IF  0$=" Y"  THEN  PRINT  #2,W$: POWER 
2950  PRINT 

2960  REM  ********  CONTINUE  ITERATIONS  ********** 

2970  PRINT  "CONTINUE  ITERATIONS  (Y/N)"; 

2980  INPUT  C$ 

2990  IF  C$="Y"  THEN  1330 
3000  GOTO  3400 

3010  REM  ********  NEW  VOLTGAES  ********** 

3020  FOR  1=1  TO  NF 

3030  VR( I )=FR( I )-C*G( 1 ,  1 ) 

3040  VI(I)=FI(I)-C*G(I,2) 

3050  NEXT  I 

3060  REM  ********  ACHIEVED  GAIN  ********** 

3070  REM  -  ELECTRIC  FIELD  FOR  ACTUAL  VOLTAGES 

3080  FOR  J=1  TO  NP 
3090  ER ( J ) =0 
3100  E I ( J ) =0 
3110  FOR  1=1  TO  NF 

3120  ER(J)=ER( J)+VR( I ) *AR( I, J)-VI ( I )*AI ( I, J) 

3130  El ( J)=EI ( J ) +VR ( I )*AI ( I , J)+VI ( I ) * AR ( I , J ) 

3140  NEXT  I 
3150  NEXT  J 

3160  REM  -  INPUT  POWER  CALCULATION  FOR  ACTUAL  VOLTAGES 

3170  POWER=0 
3180  FOR  1=1  TO  NF 
3190  P1=0 
3200  P2=0 

3210  FOR  J=1  TO  NF 

3  2  20  PI  =  P 1 +YR ( I , J ) *VR ( J ) -YI ( I , J)  *VI  ( J) 

3230  P2=P2+YI( I , J ) * VR ( J ) +V I ( J ) * YR ( I,J) 

3240  NEXT  J 
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3250  POWER=POWER+VR( I ) *P1+VI ( I ) *P2 

3260  NEXT  I 

3270  POWER=POWER/2 

3280  REM  ********  LEAST  SQUARES  ERROR  ********** 

3290  S=0 

3300  FOR  1=1  TO  NP 

3310  TT(I)=ER(I)*ER(I)+EI(I)*EI(I) 

3320  REM  -  ACHIEVED  GAIN  (DB)  FOR  ACTUAL  VOLTAGES 

3330  A( I ) =L*LOG ( GC*TT ( I ) /POWER) +LR 

3340  REM  -  DIFFERENCE  IN  ACHIEVED  AND  REQUIRED  GAIN 

3350  T(I)=R(I) -A{ I ) 

3360  REM  -  LEAST  SQUARES  ERROR  WITH  WEIGHTING 

3370  S=S+W ( I ) *T ( I ) *T ( I ) 

3380  NEXT  I 
3390  RETURN 
3400  END 
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